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A very simple formula is presented that relates the logarithm of the half-life, corrected by the 
centrifugal barrier, with the Coulomb parameter in proton decay processes. The corresponding 
experimental data lie on two straight lines which appear as a result of a sudden change in the 
nuclear shape marking two regions of deformation independently of the angular momentum of the 
outgoing proton. This feature provides a powerful tool to assign experimentally quantum numbers 
in proton emitters. 
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Nearly a century ago one of the most challenging problems in physics was the understanding of the 
Geiger-Nuttall rule, which says that the logarithm of the half-life in alpha decay is inverse proportional 
■ to the energy of the outgoing alpha particle, i. e. its Q-value. It can be asserted that the probabilistic 
interpretation of quantum mechanics started by the explanation given by Gamow to that rule as a 
t-H ■ consequence of the quantum penetration of the alpha particle through the Coulomb barrier [lj . Nowadays 
the challenge in nuclear physics is related to rare nuclei, i. e. nuclei lying very far from the stability 
line which decay rapidly by particle emission (neutron and proton drip lines) 0] . The exploration of the 
drip-lines is one of the ambitions of the projected radioactive nuclear beam facilities as e. g. the Rare 
Isotope Accelerator (RIA) 3] . 

By plotting the logarithm of the half-lives as a function of Q p , i. e. the kinetic energy of the emitted 
proton, one does not obtain a clear graphical pattern of the experimental data, as for the Geiger-Nuttall 
rule. As seen in Figure l.a 0, such a plot does not reveal any obvious trend. The reason of this disorder 
is that not only Q p but also the height and width of the Coulomb and centrifugal barriers, in which 
w ■ the proton is trapped before decaying, determine the decay width. The role of the centrifugal barrier 
__L ' in proton emission is more important than in a-decay due to the smaller proton reduced mass and also 
because in most cases the proton carries a non-vanishing angular momentum. The half-life corresponding 
to the wave with given quantum numbers can be written as [j| 
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where is the Coulomb-Hankel spherical wave, x = 2{Z — l)e 2 / (hv) is the Coulomb parameter (which 
determines the Coulomb barrier), Z is the charge number of the mother nucleus, v = hk/fi = y // 2Q p / fi 
is the velocity of the outgoing proton, a is the reduced mass of the proton-daughter system, p — kR and 
the spectroscopic function is given by [6j 

Slj (R,[3) = Y, K W?( R ^) s v]>fi'i'{ R iP) > (2) 
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where sj^ is the so-called spectroscopic factor (the particle amplitude uij within the super fluid model) 
and fij are the components of the internal wave function (the standard Nilsson function within the 
adiabatic approach). We mention here that |sy | 2 ~ 0.5. Moreover for spherical emitters it was evidenced 
a clear dependence of the spectroscopic factor upon the shell occupancy 7| . The asymptotic propagator 
matrix is given by 

K w , r (R, (3) = Hl +) ( X , p) \n {+) {R, 0)1 ^ = 8 w 5 r y + &K m , r (R, p) . (3) 

Here H^\, ^ (R, (3) is the matrix of solutions with an outgoing boundary behavior, i.e. 
'Hij^j/jt (R, P) ^b.^oc Su>5jj>Hj(Xi P)- Thus, the relation has the same form as in the spherical 
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case, but for deformed systems the spectroscopic function is a superposition of the different channel 
components /y. 

The ratio H< /sy, entering Eq. (which is just the inverse of the scattering amplitude in the 
(l,j) channel) does not depend upon the radius and this is an important test of accuracy in any coupled 
channels scheme. A good approximation of the function for physical situations is given by the WKB 
value, i. e. p| 

H l + \x,p) ~ Gi(x-.p) ~ Ci(x,p)(ctg a) 1/2 exp[x(a - sin a cos a)} , (4) 

where Gi is the irregular Coulomb function. The influence of the centrifugal barrier is fully contained in 
the function Ci, which is given by 
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V C (R) is the Coulomb potential at distance R. We will choose this distance as the matching radius, for 
which we will adopt the standard form, i. e. R — 1.2(AlJ 3 + Ap ), where Ad is the mass number of 
the daughter nucleus and A p = 1. The dependence of the centrifugal factor upon the distance, entering 
through a, is very weak around that value of R. Defining a reduced half-life as 

rp _ ZiV2 _ F( X ,P) , . 

lred ~ Cf - \ Slj {R,P)\* ' [b) 

where 

F(Xi p) — c ^9 a ex P[^x( a ~ s * n a cos a )] > (?) 
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one sees that T re d should not depend upon the angular momentum / if plotted against the dimensionless 
Coulomb parameter x- Moreover, logio T re< i is a linear function of x independently of the value of I, 
provided the velocity v, the function sy and the parameter a, which depend upon the Q-value, has 
a smooth behavior in this logarithmic scale. Notice that the value of T re d is model independent since 
the only theoretical quantity entering in its definition is the function Cf, which is a by-product of the 
barrier penetration process. But it is important to stress that this entire analysis is based upon the 
assumption that the proper value of I, that is the one that determines the experimental half- life Tin, is 
used. Otherwise, and since Ti/ 2 is strongly ^-dependent, that straight line pattern would be completely 
spoiled. 

To check the rather straightforward conclusions reached above we evaluated logio T re d in cases where 
experimental data are available, as shown in Table 1. We considered proton emitters with Z > 50 and 
angular momentum as given in Ref. 0. The orbital and total angular momenta of the proton in the 
mother nucleus labeling a given partial wave will be denoted by (Imijm)- The corresponding outgoing 
values (at large distances) will be Angular momentum conservation requires that J m = Jd + j, 

where J m (Jd) is the angular momentum of the mother (daughter) nucleus. In the case of odd decaying 
nuclei to the ground state with Jd = 0, it follows that J m = j m . If, in addition, the nucleus is spherical 
then j m = j and l m = I, where l m is the orbital angular momentum of the single quasi-proton state in the 
mother nucleus. The value of j m in Table 1 is the angular momentum of the proton outside the nucleus. 
It is often extracted from theoretical predictions. 

With the values of I, x an d logio T re d of Table 1 we produced the plot shown in Figure l.b. Amazingly 
enough, the points lie all approximately along two straight lines. The nuclei on the upper line correspond 
to the first twelve cases in Table 1, i. e. those with Z < 68. It is important to point out that the angular 
momenta in all cases of Figure l.b are as reported in the literature, except the nucleus 109 I for which 
we changed the angular momentum I = of Ref. Q to I — 2. We did this because we performed a 
coupled channels calculation showing a larger component with 1=2 (0.37) in comparison with 1=0 (0.02) 
for the 1/2+ deformed state at the Fermi level. Note that this assignment is also in agreement with the 
systematics in heavier isotopes, like e.g. in m I 8]. 

Using a fitting procedure we found that the experimental half-lives can be reproduced by the formula 

l °9w T^ ] d = a k ( X -20) + b k , 

ai = 1.31 , bi = -2.44 , Z < 68 

a 2 = 1.25 , b 2 = -4.71 , Z > 68 . (8) 
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where k — 1 corresponds to the upper line in Figure l.b. The standard errors are <j\ — 0.26 and o-i = 0.23, 
corresponding to a mean factor less than two. 

These two straight lines may have been induced either by an abrupt change in the Q-values or in the 
structure of the different emitters, or by both. The Q- value dependence affects only the function F(x, p). 
We therefore plotted this function against x m Figure 2. One sees in this Figure the same pattern as 
in Figure l.b, namely there appear two lines, but now the upper line includes only six emitters with 
Z < 68. In fact the deviations from the two lines in Figure 2 are much larger here than in Figure l.b, 
corresponding to an error of about one order of magnitude. This forces us to conclude that, although 
the pattern shown in Figure 2 corresponding to sij(/3) = const is similar to the dependence between p 
and Xi this effect can not be the only reason behind the two straight lines of Figure l.b. The only other 
source that may contribute to that alignment is an abrupt change in the nuclear structure of the emitters. 
We therefore correlated the deformation values calculated in Ref. with the determined half-lives. In 
Figure 3. a we show the deformation parameter [3 as a function of Z. At the proton drip line between Z—67 
and Z=69 occurs a pronounced change, from a large prolate shape with w 0.3 to an oblate shape with 
f3 w -0.2. These shapes are substantiated by measurements of moments of inertia. In order to analyze the 
dependence of the half-life upon the internal structure of the nucleus we plot in Figure 3.b the quantity 
logiosjj 2 as a function of the charge number Z. It is worthwhile to emphasize once again that this is a 
model independent function. 

One remarks a striking correlation between Figures 3. a and 3.b. After the jump occurring at Z=68, 
where an abrupt shape change occurs, the deformation of the nuclei lying on the lower line smoothly 
increases. One can then assert that these lines reflect two regions of nuclei separated by a sharp transition 
between the prolate and oblate regimes. 

There are two points in Figure l.b that deviate conspicuously from the upper straight line. They 
correspond to 140 Ho (1=3) and 141 Ho* (1=0). In Figure 3. a one sees that these nuclei are situated at the 
border between the two regions of deformation. 

In conclusion we have presented in this paper a simple formula for proton decay (Eq. ©) similar to the 
Geiger-Nuttall rule. This formula enables the precise assignment of spin and parity for proton decaying 
states. The only quantities that arc needed arc the half-life of the mother nucleus and the proton Q- value. 
As a function of these quantities, corrected by the centrifugal barrier (Eq. ©), the experimental data 
of proton emitters with Z > 50 lie along two straight lines. Taking into account that 140 Ho 141 Ho* are 
the only emitters that deviates from the systematics, one concludes the remarkable fact: proton emitters 
with Z > 50 follow simple systematics, which are obviously correlated with two regions of quadrupole 
deformation, as can be seen by plotting logios~ 2 versus (3. 
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Table 1 

Data used in the calculation. The spin of the proton moving inside the mother nucleus (jm), the proton 

Q-value (Q p ) and the emitter half-life (T exp ) are from Refs. J2. 13. [7^/ . The quadrupole deformation 
parameters [3 are from Ref. and the quantities I, \, a, T rec i and Sij are as explained in the text. All 
half-lives are in seconds. Stars in the emitters indicate excited states. 
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FIG. 1: (a) Logarithm of the experimental half-lives corresponding to proton decay as a function of the Q-value. 
The data are taken from Ref. Q. 

(b) Values of logio T re d, Eq. @, as a function of the Coulomb parameter \i The numbers labelling the different 
symbols correspond to the I- values of the outgoing proton. The two lines are computed according to Eq. JHJ . 
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FIG. 2: The logaritm of F(x,p) (s/fm), Eq. (JJJ, as a function of \. The two straight lines are average values to 
guide the eye. Filled symbols correspond to the first twelve cases of Table 1. 
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FIG. 3: (a) Dependence of the quadrupole deformation parameter /3, as given in Table 1, upon the charge number 
Z. Filled symbols correspond to the first twelve cases of Table 1. Notice that (3 is the same in this scale for the 
isotopes Cs and also for the two states in Ho. 

(b) The logarithm of the spectroscopic function s ; ~ 2 , defined by Eq. Q and given in the last column of the Table 
1, versus the charge number Z. 



